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Abstract 

Quasipatterns have been found in dissipative systems ranging from Faraday waves 
in vertically vibrated fluid layers to nonlinear optics. We describe the dynamics of 
octagonal, decagonal and dodecagonal quasipatterns by means of coupled Ginzburg- 
Landau equations and study their stability to sideband perturbations analytically 
using long-wave equations as well as by direct numerical simulation. Of particu- 
lar interest is the influence of the phason modes, which are associated with the 
quasiperiodicity, on the stability of the patterns. In the dodecagonal case, in con- 
trast to the octagonal and the decagonal case, the phase modes and the phason 
modes decouple and there are parameter regimes in which the quasipattern first 
becomes unstable with respect to phason modes rather than phase modes. We also 
discuss the different types of defects that can arise in each kind of quasipattern as 
well as their dynamics and interactions. Particularly interesting is the decagonal 
quasipattern, which allows two different types of defects. Their mutual interaction 
can be extremely weak even at small distances. 



1 Introduction 



Since the discovery of quasicrystals in 1984 [1], much attention has been paid 
to the properties of these materials. In contrast to perfect crystals they lack 
periodicity, but preserve long-range orientational order. Due to this lack of 
periodicity quasicrystals may have non-crystallographic rotational symmetry 
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and, in fact, materials have been found with five, eight, ten or twelve- fold 
rotational axes. 

In dissipative systems the possibility of quasiperiodic structures, or quasi- 
patterns, was also suggested some time ago [2]. Since then they have been 
observed in Faraday waves with one- and two-frequency forcing [3,4] and in 
nonlinear optics [5]. Marangoni convection with a deformable interface has 
also been suggested to support quasipatterns [6] , although they have not been 
observed so far. 

In Faraday waves a great variety of patterns has been observed. They include 
superlattices, rhombic states, oscillons, as well as quasipatterns [7-9]. In or- 
der for superlattices and quasipatterns to be stable the mutual suppression of 
plane-wave modes of different orientation has to be sufficiently weak. This can 
be due to a few, somewhat different mechanisms. When forced with two fre- 
quencies, the system can exhibit simultaneously instabilities at two different 
wavelengths and the quasipatterns and superlattices appear near this bicritical 
point. The angle between the wavevectors of the destabilizing modes depends 
on the ratio of the two wavelengths. It determines whether the resulting pat- 
terns corresponds to a superlattice or a quasipattern. A mechanism involving 
two length scales was suggested some time ago for quasicrystals [10] and stud- 
ied in dissipative systems by means of a modified Swift-Hohenberg equation 
with two marginal modes [11]. The second length scale need not be associated 
with a proper instability. It can be sufficient that one of the modes is weakly 
damped. The interaction of the unstable mode with this damped mode can 
then lead to a reduction in the mutual suppression of modes subtending a 
certain angle [12]. Alternatively, the interaction with the damped mode can 
strongly enhance the saturating self-coupling term which then effectively leads 
to a reduction of the competition over a wide range of angles [13,14]. This fa- 
vors patterns with a large number of modes. Calculations of the coefficients 
of the amplitude equations for Faraday waves suggest that this is indeed the 
case [15,16]. 

In nonlinear optics the optical field in a nonlinear cavity can be rotated in 
order to obtain a structure with the desired number of modes [5]. In these 
systems it is also possible to arrange patterns with two different spatial scales, 
leading to complex quasicrystaline structures [17-19]. 

The relative stability of perfect quasipatterns has already been addressed [2]. 
In the present paper we are interested in the stability of n-fold quasipatterns to 
sideband perturbations, in particular to slowly varying modulations. We will 
assume that the physical fields can be expanded as a sum of Fourier modes 
rotated by 1-njn relative to each other (similar to the density wave picture 
of quasicrystals), with slowly varying amplitudes in space and time. Using 
symmetry arguments we determine a set of coupled Ginzburg-Landau equa- 
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tions for these amplitudes. Due to the lack of periodicity, a rigorous deriva- 
tion from the basic equations using a center manifold reduction has not been 
possible so far. However, we will consider the Ginzburg-Landau equations as 
phenomenological models and aim to extract conclusions that can be verified 
experimentally. 

As is well known in the study of quasicrystals, the long-wave dynamics of 
the system is governed by two types of modes: phonons and phasons. The 
former ones are marginal modes arising from spatial translation symmetries. 
The phasons, on the other hand, are additional marginal modes characteris- 
tic of quasicrystals and appear because of the quasiperiodic nature of these 
structures. As we will see later, in octagonal and dodecagonal quasipatterns 
they have a very simple geometrical interpretation, since these quasipatterns 
can be viewed as two superimposed square (hexagonal) lattices rotated by 
27r/8 (27r/12) . The phason modes correspond to relative translations between 
the two lattices. It is worth mentioning that when the angle between the two 
lattices satisfies certain conditions, the whole structure becomes periodic with 
a wavelength larger than that of the individual lattices: this is known 
perlattice [20]. In this case, extra resonance conditions among the modes are 
satisfied [21], and the phason modes become damped and can, in principle, be 
eliminated in a long-wave analysis. 

The paper is organized as follows. In the next three sections we study octag- 
onal, decagonal and dodecagonal quasipatterns arising from steady bifurca- 
tions. In each case we start with the amplitude equations and study various 
simple solutions and their relative stability. Our main focus is the long-wave 
dynamics of the quasiperiodic solutions. We derive coupled equations for the 
phonon- and phason-type modes, obtaining the diffusive counterpart to the 
elastic equations for a quasicrystal [22]. We then calculate the long- wave sta- 
bility limits for sideband perturbations and investigate numerically the be- 
havior arising from the instabilities. In section 5 we discuss the different types 
of defects that appear in each kind of quasipattern and their dynamics and 
interactions. Conclusions are presented in section 6. 



2 Octagonal Quasipatterns 

2.1 Ginzburg-Landau Equations 

We consider a quasipattern composed of four modes, ip = rj J2 n =i Ai( x , t)e lk "' x + 
c.c. + h.o.t, where 77 <C 1 and the amplitudes Aj(x, t) satisfy the equations (af- 
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Fig. 1. Fourier modes composing the a) octagonal, b) decagonal and c) dodecagonal 
quasipatterns. Solid vectors denote the modes described by Eqs. (1), (39), and (85, 
86), respectively. 

ter rescaling): 



d t A. = Mi + (n* • V) 2 Aj - A,i \A,\ 2 + u\A l+2 \ 2 + j(\A i+1 \ 2 + \A l+3 \ 2 ) , (1) 



with \x the control parameter and the unit vector in the direction of the 
wavenumber kj (see Fig. la). The amplitudes depend on slow scales x = rjx 
and t = rft. Here and in the following the indices are repeated cyclicly with 
period 4. Thus, A 5 corresponds to Ai, etc. The coefficients v and 7 measure the 
interaction between modes that subtend an angle of 7r/2 and 7r/4, respectively. 

Eqs. (1) have gradient structure 

d t Ai = -U, (2) 
dAi V ; 



and can be derived from the Lyapunov functional 



4 

T = J J dxdyF = J J dxdy ^ 

i=l 



-fi\Ai 



(n,-V) 2 A| 2 + -|A 



+-f\Ai\ 2 \A i+1 \ 2 + -\Ai\ 2 \A 



i+2\ 



(3) 



Therefore, the dynamics of the system are relaxational. 



In order to calculate the different solutions we write Aj = cije 1 ^ , with a,j real. 
Splitting into real and imaginary parts yields 



dttti = /Mi - di[a 2 + ua 2 i+2 + 7(a 2 +1 + a 2 +3 )], (4) 

d t <j) t = 0. (5) 

The amplitudes can be considered to be real (0i = • • • = 04 = 0). There are 
five different kinds of simple solutions [2] (see Fig. 2): 
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Fig. 2. Octagonal case. Different types of solutions of Eqs. (1), obtained by recon- 
structing the physical field tp from the amplitudes A^. a) Rectangles, b) squares, c) 
ID quasipattern and d) octagonal quasipattern. 



(1) Rolls: a 2 = a 3 = a 4 = and oq = ^JJL. The value of the Lyapunov func- 
tional for rolls is: F Ro u = — /i 2 /2. 



(2) Squares: a 2 = a 4 = and a x = a 3 = \J fi/ \v + 1), F Sq = —fi 2 /( y u+ 1). 

(3) Rectangles: a 3 = a 4 = and a x = a 2 = y / /Vv7 + 1), F Rect = -/i 2 / (7+ !)■ 

(4) One-dimensional quasipattern: a 4 = and 



a 2 



///(l + t/-27) 
1 + v - 2 7 2 : 



ai = a 3 



; Mi -7) 

1 + v - 2 7 2 



(6) 



1DQ 



fx 2 3 



47 



2 1 + v - 2 7 2 



(7) 



It is quasiperiodic in the direction of k 2 and periodic along k 4 . As is typ- 
ical for solutions with submaximal isotropy, this solution exists only over 
a limited range of the nonlinear coefficients, (1 + v — 27) (1 — 7) > 0. It 
is always unstable [2]. 
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(5) Quasipattern: a± — • • • — a& — R, with R satisfying, 

R= \jl + u + 2 1 Fq = ~1 + v + 2 1 (8) 

The conditions for linear stability of these solutions are given by 

(1) Rolls: v > 1, 7 > I. 

(2) Squares: 7 > (1 + u)/2, -1 < v < I. 

(3) Rectangles: v > 1, — 1 < 7 < 1. 

(4) Quasipattern: -(1 + v)j2 < 7 < (1 + u)/2, -1 < v < I. 
2.2 Longwave analysis 

We assume that the conditions for the two-dimensional quasipattern to be 
stable with respect to homogeneous perturbations are satisfied and study its 
stability to sideband perturbations. In particular, we study the stability of 
the pattern as a function of its wavenumber. A perfect quasipattern with 
wavenumber k slightly different from critical is given by Ai = i?e lljni ' x , with 
R=\A 1 \ = \A 2 \ = \A 3 \ = \A 4 \ = v /( / u-g 2 )/(l + z/ + 2 7 ) and rjq = |k - k c |. 
We expand around this solution, A { = (R + r^e 1 ^'^^. Considering first 
only homogeneous perturbations (r, = Tj(t), fa = fa(t)) and separating real 
and imaginary parts, the linearized equations for the perturbations are: 

d tri = -2R\ - 2uR 2 r i+2 - 2 1 R 2 {r l+l + r l+3 ), (9) 
d t fa = 0. (10) 

For the amplitude perturbations there are three eigenvalues, two corresponding 
to one-dimensional eigenspaces, 07^ = — 2i? 2 (l + z/ + 27), for = [r± = l,r 2 = 
1,7-3 = l,r 4 = l]/2, and cr? 2 = -2R 2 (1 + v - 2 7 ) for v? 2 = [-1, 1, -1, l]/2. 
The other eigenvalue <Jf s ,r 4 = —2R 2 (1 — v) corresponds to a two dimensional 
eigenspace spanned by Vf 3 = [—1, —1, 1, l]/2, v? 4 = [—1, 1, 1, — 1]/2. 

The four phases are marginal modes. Physically it is more convenient to re- 
arrange them into phonon and phason modes. The first two are related to 
translations in space and can be written as 

uJ B = [0i = 1, 02 = cos(2tt/8), 03 = cos(4vr/8), 4 = cos(6tt/8)], (11) 
= [0, sin(27r/8), sin(47r/8), sin(67r/8)]. (12) 
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There still remains a two-dimensional subspace orthogonal to this one. An 
orthonormal basis for it is given by 



u T = [1, - cos(2tt/8), cos(4vr/8), - cos(6tt/8)], (13) 



<pi 

T 



u 1 ^ = [0, - sin(27r/8), sin(47r/8), - sin(67r/8)]. (14) 



This choice of the phason modes corresponds to relative translations (in the x- 
and y-directions, respectively) between the two square lattices that compose 
the quasipattern. It is worth mentioning that the phason modes do not trans- 
form as a vector. In fact, under rotations by an angle 9 the transformation 
of the phason field <p = (ipi, <p 2 ) corresponds to that of a regular vector for a 
rotation by an angle 59 (cf. (30.31)). 

When spatial modulations are included the perturbation equations become 



d tri = -2qR{h t • V)0, + (n, • V) 2 r, (15) 
-2R\ - 2uR 2 r i+2 - 2 1 R 2 (r l+l + r l+3 ), 

^ = (n l -V) 2 0, + ^(n i -V)r,. (16) 

Now the phase modes are no longer marginal, but exhibit diffusive dynamics. 
In order to study these long-wave dynamics we define a small parameter e and 
introduce slow time and space scales: T = et, X = e 1 / 2 x. The amplitude and 
phase perturbations are expanded as: 



r^XT) = e£r}(X,7K ( 17 ) 
&(X, T) = e 1/2 [0,(X, T)u\ x + (f> y (X, T)ui y , 

+^(x,T) lt ; i + P2 (x,T)«;j. (is) 



Inserting the expansion into Eqs. (15), (16), we obtain at order e 1 / 2 the eigen- 
values for the homogeneous perturbations, at order e a relation between the 
stable and marginal modes of the type fj = fi(V(f> x , V0 y , V</?i, V</?2), and at 
order e 3 / 2 a solvability condition for the marginal modes. This gives us the 
slow, longwave dynamics. In component form the final equations are: 

d T <j )x = D 1 V 2 ( j )x + (D 2 - DjdxW ■ $) (19) 

+D 3 (d 2 x ip 1 - 2d 2 XY ip 2 - dyifii), 
d t <p y = D^ 2 ^ + (D 2 - D^dyiV ■ 0) (20) 

+D 3 (-d 2 x ip 2 - 2d 2 XYVl + d 2 Y <p 2 ), 
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d m = - \{D A + D 5 )VVi + {Da - D 5 )(d 2 x <p 2 - d Y <p 2 - 2d 2 xyVl ) 



+D 3 (d x (fr x - 2d XY (f)y - d Y (f) x ), 
1 r, 
2 



= - [(D 4 + D 5 )V 2 ip 2 + {D 4 - D 5 )(d 2 xVl - <9 2 ^ + 2d 2 XY <p 2 )] 
+Dz(-d 2 x <P y -2d 2 XY <p x + d 2 Y <P y ), 

where the values of the coefficients are given by 



(21) 
(22) 



1 a 2 

D i = j-—, (23) 
4 u 

D 2 = \- q -- 2 -f=D 1 + D> 2) (24) 
4 u vi 

D 3 = D 1 , (25) 
D A = D u (26) 

D 5 = 3 -- q -- 2 -^ = D 1 + D>, (27) 
4 u v 2 

with uj = 2R 2 (l + u + 2 7 ), w 2 = 2R 2 (1 + v - 2 7 ), u = 2R 2 (l - v), and R 
given by Eq. (8) with \i replaced by ji — q 2 . Eqs. (19)- (22) are the diffusive 
analogue of the elastic equations for octagonal quasicrystals, and their form 
can be deduced directly by means of symmetry arguments. 

A more compact notation, in which the symmetries become more evident, can 
be achieved by writing the phonon and phason modes as complex fields. Let 
= <Px + i0 y , V 9 = Vi + ' 1( P2 an d V = dx + idy. Then the former expressions 
become: 



d t (f) = £>i|V| 2 + (D 2 - Di)V(V0 + V0) + D 3 e ia W 2 ip, (28) 
dtip = \{D 4 + £) 5 )|V|V + ^(A - A,)V 2 y5+ /V Q V 2 (29) 

The angle a depends on the basis that is chosen for the phonon and phason 
modes. Since the specific choice does not have any physical relevance, a does 
not appear in the dispersion relation. In our case, a = 0. Now it is easy to 
see that these equations are invariant under rotations of nn/4. In fact, under 
such a rotation — > e n7n / 4 0, V — > e n7rl / 4 V, (f — > e 5ri7rl / 4 (p, and consequently 

V 2 -> e- 2ri " i/4 vV n7ri/4 = e- 3n7ri/4 V 2 0, (30) 
VV -> e - 2 ™ /4 vV 5Mi/ V = e rai/4 vV. (31) 



Thus, while V transforms as ip, V y? transforms as 0. 
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Fig. 3. Stability diagrams for a) v = 7 = 0.5, b) v = 0.7, 7 = 0.4. The dotted, 
dashed, and dot-dashed lines correspond to D' 2 = 0, D' 5 = and D\ = 0, respec- 
tively. The circles are obtained solving the full dispersion relation associated with 
Eqs. (15), (16). 
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In order to calculate the stability boundaries, normal modes must be consid- 
ered: = 0oe iQ ' x , (p — y?oe iQ ' x . From Eqs. (19)-(22), it can be shown that the 
stability curves are given by the expression: 



2D X D 2 D A D^ = D 2 ^ + D 2 )(D 4 + D 5 ) - 2D 2 3 ], 



(32) 



or, using Eqs. (23)-(27), 

D\D' 2 D' b = 0. (33) 



Thus, two eigenvalues become zero when Di — and the other two when 
D' 2 = and D' 5 = 0, respectively. From this, the stability curves are given by: 



D x = 



3 + Z/ + 47 2 

»i = — ; 1 , 



D' 2 
D> 








3g 2 , 



3 + 3z/ + 27 2 
1 + v — 27 



(34) 
(35) 
(36) 



It is interesting to note that the instability corresponding to D' 2 = is given 
by the usual value for the Eckhaus curve. Typical stability diagrams are shown 
in Fig. 3. When v — 7 the problem is degenerate and /ii = /i 5 (Fig. 3a). Which 
instability occurs first depends on the values of the coefficients v and 7. For 
7 > 0, v < 7; v + 7 < 0, 7 < 0; and z/ + 7>0, z/>7 the first instability is 
that corresponding to D' 5 , D' 2 , and D ± , respectively. The third case is shown 
in Fig. 3b. The symbols denote the results of a stability analysis using the full 
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dispersion relation without the long-wave approximation. They show that the 
relevant instabilities are long-wave. 

In general, the eigenvalues depend on the angle of the perturbation. When 
9 = &Yctan(Q y /Q x ) = rnr/i there is an eigenvalue that is always marginal. 
(From Eqs. (15), (16) it is easy to see that the growth rate of a perturbation 
to the phase fa in a direction perpendicular to hi, is zero). This means that 
the perturbations perpendicular to each set of rolls evolve on a still slower 
time scale that is not captured with Eqs. (1). This situation is equivalent to 
that of the zig-zag instability of rolls or squares [23,24]. In order to resolve this 
degeneracy we must take Newell- Whitehead-Segel derivative terms in Eqs. (1), 
using the replacement 

(n.-V) 2 -> [(n, ■ V) - WV 2 ] 2 , (37) 

where 5 is a small coefficient, whose size depends on the distance from thresh- 
old, 5 = 0(r]). 

Assuming that the phase of each of the modes only depends on the direction 
perpendicular to that mode, it is easy to see that the equations decouple. 
Introducing a super-slow time, dr = e 4 <9r 4 , and taking q to be small, q = e 2 q, 
the usual nonlinear phase equation for the zig-zag instability is obtained: 

d Ti fa = 2q5(n ■ Vffa - 5 2 (f, • V)% + 65 2 [(f, • V)fa]\f t ■ V) 2 fa, (38) 

where fj is the unit vector perpendicular to fij. When q < the pattern is 
unstable. As the coefficient in front of the cubic term is positive, the insta- 
bility is supercritical. In order to confirm this we have performed numerical 
simulations of Eqs. (1) using a fourth-order Runge-Kutta method with an in- 
tegrating factor that computes the linear derivative terms exactly in Fourier 
space. In Fig. 4 the evolution of this zig-zag instability is shown, starting 
with a perturbation of the quasipattern in the direction perpendicular to h^. 
A l = Re iqx (l + 0.1icos(47n//L)), Aj = Re iq&i ' x , j ^ 1. As expected, the 
instability saturates, resulting in a distorted quasipattern. 



3 Decagonal Quasipatterns 

3.1 Ginzburg- Landau Equations 

For the decagonal quasipattern we consider the expansion (cf. Fig. lb) ip = 
vT,n=i A n (x,t)e t]ln ' x + c.c. + h.o.t.. The amplitudes Aj(x, t) = rje'** now satisfy 
the equations 
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Fig. 4. Reconstruction of the octagonal quasipattern from the amplitudes Ai for the 
zigzag instability at the times a) t=0, b) t=1000, and c) t=3000. The simulations 
are done with 64 x 64 Fourier modes and the values of the coefficients: fi = 5, 
v = 0.7, 7 = 0.4, q = -0.05, L = 25, k c = Wk min and 5 = 0.1. 



d t A = A*A + (n* • V) 2 A - -4i[|A| 2 + KIA+i| 2 + |A+4| 2 ) 



+7(1^4 



j+2 



'A+3| 2 )] + aAi +1 A i+2 A i+3 A i+4 , 



+ 



(39) 



where // is the control parameter and v and 7 measure the interaction between 
modes subtending an angle of 27r/5 and 47r/5, respectively. The indices are re- 
peated cyclically with period 5. Although the quartic term is higher order than 
the others, we have included it in Eq. (39) to account for the resonant interac- 
tion among the five wavevectors, Yjj=\ k? — (see Fig. lb). Without this term 
there would be a spurious one-parameter class of solutions, parameterized by 
the global phase, $ = Y^=i <j>ji an with the same energy. The resonance term 
lifts this degeneracy by rendering the global phase (slightly) damped. 



The Lyapunov functional for the dodecagonal case is given by 



dxdy 



i=i 



-^A^ + ^-VfA^ + ^A 



2 A 12 , ,4 hz/ |4.|2M. 12 



H+l I 



+7|Ai| 2 |A i+2 | 2 ] - a(A 1 A 2 A 3 A 4 A 5 + c.c). 



(40) 



Writing again Aj = we now have: 



dttti = fiai - a^a] + z/(a 2 +1 + a 2 +4 ) + 7 (a 2 +2 + a 2 +3 )] 

+aa i+ ia i+2 ai +3 a i+ 4COs(<I>), (41) 
aid t <pi = -aaj + iai +2 ai + 3ai + 4 sin($), (42) 

with $ = X^=i 0j the global phase. For a quasiperiodic solution ai — ■ • • — 
a 5 = R and Eq. (42) becomes: 

d t $ = -5a J R 3 sin($). (43) 
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Fig. 5. Decagonal Case. Different types of solutions of Eqs. (39). a) Rolls, b) rectan- 
gles Ri, c) rectangles R2, d) and e) ID quasipatterns H± and H2, and f) decagonal 
quasipattern. 

There are two solutions, one stable ($ = 0), the other unstable ($ = n). We 
will consider the former one. 



Because the resonant term involving a is quartic, it cannot appear in sys- 
tems with the reflection symmetry Aj — > — which arises, for instance, in 
Boussinesq convection or in Faraday waves with one-frequency forcing. In 
that case (a — 0) $ is arbitrary to the order considered in Eq. (39) and, for 
$ 7^ 0, 7r, solutions with pentagonal rather than decagonal symmetry can bi- 
furcate from the trivial state. The situation is analogous to that of the triangle 
solutions in the case of hexagonal symmetry [25] . The global phase gets fixed 
by the higher-order resonance term Aj|Aj + i| 2 |Aj + 2| 2 |Aj + 3| 2 |Aj + 4| 2 , which leads 
to <9j$ ~ sin(2$), implying $ = mi/2 for the regular pentagonal solution. We 
will not consider these solutions in the following and take 0i = • • • = 5 = 0. 

There are six different kinds of simple solutions (see Fig. 5): rolls, two types 
of rectangles, two types of one-dimensional quasipatterns and the decagonal 
quasipattern. They are given by: 



(1) Rolls: a 2 = • • • = a 5 

(2) Rectangles (-Ri): cq 

-/iV^ + i). 

(3) Rectangles (-R2): o 2 



= and cq = y/Jl, F R = 
= a 2 = 05 = and 03 

= a 3 = a 5 = and cq 



-/i 2 /2. 

V^+l), F Rl = 
VV(7 + 1), F R2 = 



= CI4 = 



= a 4 = 



12 



-^7(7 + 1)- 

(4) One- dimensional quasipatterns (Hi): a 2 = a 5 = and 



/i(l + z/-2 7 ) A*(l -7) 

01 = V l + »/-2f» ' ° 3 = ° 4 = V l + ^-2T" ^ 



^ 1/ + 3 - 47 



2 1/ - 2 7 2 + r 

They exist provided (1 — 7 )(1 + v — 2 7 ) > 0. 
(5) One-dimensional quasipatterns (H 2 ): a 3 = a 4 = and 



(45) 



0l = V^ a2 = a5 = VT^ 0~~2' 46 

V 1 + 7 — 2z/ 2 V 1 + 7 ~~ 2i/^ 



/i J 3 + 7 - 4z/ 
^ 2 ~ ~ 2 1 + 7 - 2z/ 2 



(47) 



They exist provided (1 — + 7 — 2z/) > 0. Both, Hi and are 
quasiperiodic in one dimension and periodic in the other. 

(6) Quasipattern: ai — • • • a 5 — R, with R satisfying 

= /i - (l + 2z/ + 2 7 ) J R 2 + ai? 3 . (48) 

The analytic solution is quite complicated. For a — it simplifies to 



*=\/ttS =- F « = -5irfe- (49) 

(7) In addition, when ct = there is another solution, with 05 = and 



V I + 7 + 1/ + 71/ — i/ 2 — 7 ^ 

/ gj -7) 

a 2 = a 3 = W — ; ; ~ 2 2 - 51 

V 1 + 7 + z/ + 7 z/ — v l — '-f 2 

It is also quasiperiodic in the two spatial dimensions and exists if (1 — 
I /)(l- 7 )>0. 

The conditions for linear stability are: 

(1) Rolls: v > 1, 7 > 1. 
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(2) Rectangles Ri. -1 < v < 1, 7 > 1. 

(3) Rectangles i? 2 : z/ > 1, — 1 < 7 < 1. 



(4) ID quasipattern Hi. 



27- 1< v < 1, 

1 + V5 1-V5 n 
1 - z/ — 7 — - — < 



(52) 
(53) 



(5) ID quasipattern H 2 : 



2i/-1<7<1, (54) 
1 — z/ 7 <0 (55) 



(6) The quasiperiodic solution given by Eqs. (50), (51) is always unstable. 

(7) Decagonal quasipattern: 



1 + 1 a/5 
1 - v — 7 V Ra > 0, 



2 ' 2 

l-VE l + VE 
-7- 



1 - z/ 

2 2 
2(l + 2z/ + 2 7 ) -3ai?> 0. 



+ i?a > 0, 



(56) 

(57) 
(58) 



When (58) is violated the quasipattern undergoes a saddle-node bifurca- 
tion that generates an unstable branch with larger amplitude (see Fig. 6). 
For these amplitudes the quartic term is of the same order as the other 
terms in Eq. (39), which is inconsistent with the amplitude expansion. 
We restrict ourselves therefore to the range 3aR -C 2(1 + 2v + 27). 

For R = Eqs. (56,57) are complementary to the stability condi- 
tions (53) and (55) for the one- dimensional quasipatterns. Thus, at onset 
one of the one-dimensional quasipatterns can be stable while the two- 
dimensional quasipattern is unstable (or vice versa). For a > 0, the two- 
dimensional quasipattern can gain stability with increasing /i while the 
one-dimensional quasipattern is still stable. This bistability is expected to 
persist even if higher-order terms were included in the Ginzburg-Landau 
equations (39), since they would contribute only terms of the order R 2 to 
the stability conditions. For a small range of parameters v and 7 the sta- 
bilization of the two-dimensional quasipattern can occur for amplitudes 
for which these higher-order terms would still be negligible. 
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R 




Fig. 6. Sketch of the bifurcation diagram for the decagonal quasipattern, as given 
by Eq. (48) (solid and dotted lines representing stable and unstable branches, re- 
spectively). The saddle- node bifurcation occurs at amplitudes for which (39) is not 
valid any more. 

3.2 Longwave analysis 

Proceeding as in the case of the octagonal quasipattern we obtain for the 
linearized perturbation equations 

5 

d t n = -2 gj R(n, • V)0, + (n, • V)\ + i? 3 a(E r j ~ 2r i) (59) 

i=i 

-2R\ - 2uR 2 (r i+1 + r^) - 2 1 R 2 (r t+2 + r,_ 2 ), 

2a ~ 5 
d t <f>i = (n, • V) 2 0, + -f (n 4 • V)n - aR 3 £ 0„ (60) 

H 3=1 

with R satisfying 

= (/i - q 2 ) - (1 + 2v + 2 7 )i? 2 + aR 3 . (61) 



For the amplitude perturbations there are three eigenvalues. One eigenvalue, 
a R = R 2 (3aR — 2 — 4z/ — 47), corresponds to a one-dimensional eigenspace 
spanned by 

v T H = -=[n = l,r 2 = l,r 3 = l,r 4 = l,r 5 = 1]. (62) 
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The other two eigenvalues are given by a T 2,3 = R 2 {-2 - 2aR + ± VE) + 
7(1 =F v^o") and each corresponds to two two-dimensional eigenspaces. Four 
orthonormal vectors spanning these spaces are 



-(i + VE)/2, (i + VE)/2, -1,0,1 



; T 3 



(63) 



-(1 - V5), -(1 - VE), -(1 + v/5), 4, -(1 + V5)] , (64) 



corresponding to 0"T2 and 



T 

v T2 = - 


1 








T 

v T4 = 1 


rr 

V 40 


-(1 



-1,(1- V5)/2, -(1-^/2, 1,0 



(65) 



1 - V5), -(1 + >/5), -(1 + V5), -(1 - >/5), 4] , (66) 



to a T3 . 

The global phase is a stable mode with a one-dimensional subspace spanned 
by m$ = [</>i = 1, . . . , 05 = l]/\/5. As with the octagonal quasipattern, there 
are four marginal phase modes. The two corresponding to translations in space 
can be written as 



M ^ = V 5 ^' cos ^ 27r / 5 )' cos ( 47r / 5 )> cos(6tt/5), cos(8tt/5)] , 



M ^ = y- [0,sin(27r/5),sin(47r/5),sin(67r/5),sin(87r/5)] 



(67) 
(68) 



The two phason modes can be written as 



= 4/-[l,cos(67r/5),cos(127r/5),cos(187r/5),cos(247r/5)], (69) 



ul 2 = J- [0,sin(67r/5),sin(127r/5),sin(187r/5),sin(247r/5)] 



(70) 



Under rotations the phason mode ip = (ip±, ^2) changes with twice the rotation 
angle. 



The expansion of the perturbations includes now the global phase, 
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(x,t)=eJ2r}(e 1 ^,etH j (71) 



&(x, t) = e^^x, + y (e 1/2 x, et)< 

+y?i(e 1 / 2 x, et)< + v^ 2 (e 1/2 x, et)<] + e$ (e 1/2 x, et)*4, (72) 

and the longwave equations can be written as: 

d t 4> x = D 1 V 2 <t> x + D 2 d x (W ■ $) + D 3 (d 2 x2ipi + 2d 2 xyV2 - d 2 m ), (73) 

d t <t> y = D 1 V 2 <P y + D 2 d v (V ■ 0) + D z {d 2 x ,y l - 2d 2 xyVl - d 2 2 <p 2 ) , (74) 

= D 4 V V + D 3 (^ 2 ^ - Zd 2 xy 4> y - c^), (75) 

= D 4 V 2 <p 2 + D 3 (d 2 x ^y + - ^v)- ( 76 ) 

The values of the coefficients are given by 



(77) 

D2 = ^~^r, (78) 

(79) 

1 a 2 

D 4 = -- — - — = D± + D 4 , with D 4 = - - — , (80) 
2 ui u 2 4 -u 2 

with mi = R 2 (2(Ra + 1) - z/(l + v 7 ^) -7(1 - v 7 ^)), w 2 = R 2 (2(Ra + 1) - - 
v 7 ^) - 7(1 + V^)), w = ^ 2 (2(1 + 2z/ + 2 7 ) - 3Ra), and i2 a solution of Eq. (61). 

In complex form (0 = (p x + i(p y , ip = (fi+i(fi 2 , V = d x + id y ) the phase equations 
read 



dt<f> = A|V| 2 + ^ 2 V(V0 + V0) + D 3 e ia V% (81) 
^ = D 4 |V| V + D 3 e ia V 2 0. (82) 



Again, the angle ct depends on the relative orientation between the phonon 
and phason modes. For the basis chosen, a = 0. 

Considering normal modes = 0oe iQ ' x , ip = y?oe lQ ' x , the eigenvalues now 
become: 



1 




4 ~ 


Ml' 


1 




2 ~ 


J 


£>i, 




1 


g 2 g 2 


2 ~~ 


Ml U 2 



1 

^1,2 = "2 



D 1 + D 4 ± y J(D 1 - D 4 f + 4D 2 ] Q 2 , (83) 
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-1.0 



0.0 
-0.75 



0.0 1.0 -0.75 0.00 

Wavenumber q Wavenumber q 

Fig. 7. Stability diagrams for a = 1 and a) v = 7 = 0.7, b) v = 0.5, 7 = 0.9. The 
dash-dotted line represents the line at which the 2-d quasipattern becomes stable. 
Below this line, the 1-d quasipattern H 2 is the only stable solution. 



0.75 



03,4 = 



D x + D 2 + D 4 ± J(D 1 + D- 



D A f + AD 2 



Q 2 



(84) 



Typical stability diagrams are shown in Fig. 7. When v — 7 the eigenval- 
ues become degenerate (_D 4 = 2Di and three eigenvalues go through zero at 
the curve u\ = u 2 = 4g 2 .) In Fig. 7b the one-dimensional quasipattern per- 
sists amplitude-stable beyond the dashed-dotted lines (cf.(55)). We have not 
investigated its stability with respect to sideband perturbations. 



4 Dodecagonal Quasipatterns 



4-1 Ginzburg-Landau Equations 



Dodecagonal quasipatterns can be thought of as a combination of two ro- 
tated hexagon patterns. The amplitudes A(x, t) corresponding to the modes 
indicated in Fig. lc satisfy the equations: 

d t Ai = fiAi + (n; • VfAi + aA I+2 A i+4 - Ai[\Ai\ 2 + v{\A i+2 \ 2 + \A i+A \ 2 ) 

+ 1 (\A l _ l \ 2 + \A^\ 2 ) + 2{3\A l+l \\ 2 = 1,3,5, (85) 

dtAi, = fiAi + (n; • WfAi + aA i+2 A i+4 - A[\A\ 2 + u(\A i+2 \ 2 + \A i+4 \ 2 ) 

+ 1 (\A l+1 \ 2 + \A l+3 \ 2 ) + 2p\A l _ 1 \ 2 }, 2 = 2,4,6. (86) 

Here the indices are cyclic with period 6, /i is the control parameter and is, 7 
and (3 measure the interaction between modes subtending angles of 2n/3, tc/6 
and 7r/2, respectively. The two hexagonal sub-lattices imply two resonance 
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Fig. 8. Dodecagonal Case. Several of the solutions of Eqs. (85), (86). a) Rectangles, 
b) hexagons, and d) dodecagonal quasipattern. 

conditions + k 3 + k 5 = and k 2 + k 4 + k 6 = 0. Correspondingly, there are 
two global phases: $i = <f>i + 3 + 5 and $ 2 = 2 + 4 + 06- 

The Lyapunov functional can be written as: 



J J i=i 



6 r 1 



fj,\Ai\ 2 + \(hi- V) 2 A| + + v\Ai\ \ A i+2\ 



-a(A 1 A 3 A 5 + A 2 A 4 A 6 + c.c). (87) 



In terms of the magnitudes a, and phases 0j Eqs. (85,86) can be written as 

d t cii = n<ii + aa i+2 a i+A cos($) 

-a t [a 2 + u(a 2 +2 + a 2 +4 ) + 7 (a? ±1 + a 2 i±3 ) + 2f3\a l±1 \ 2 }, (88) 
aid t <t>i = -aa i+2 a i+A sin($), (89) 

where $ = $1 or $ = $ 2 for % odd and % even, respectively. For hexagonal or 
dodecagonal patterns in which the non-zero amplitudes are fixed, dj = R, Eq. 
(89) becomes 

$$ 1>2 = -3a-Rsm($i )2 ). (90) 



Among the four solutions $i 2 = and <3>i j2 = 7r only the one with $! = $ 2 = 
is stable. In the following we assume the amplitudes Aj to be real and take 
0i = ■ • • = 6 = 0. Now the simple solutions include rolls, rectangles, squares, 
hexagons, mixed modes, one-dimensional quasipatterns and the dodecagonal 
quasipattern (see Fig. 8). They are given by: 

(1) Rolls: a 2 = ■ ■ ■ = a 6 = and a x = ^//Z, F R = —fi 2 /2. 
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(2) Rectangles (R): a 2 = ■ ■ ■ = a 5 = and ai = a 6 = y / /V(7 + 1)> -^Ri = 

-^7(7-i). ' 



(3) Squares (S*): a 3 = • • • = a 6 = and a x = a 2 = y / V*/(2/3+ 1), i* 1 ^ 
-// 2 /(2/3 + l). 

(4) Hexagons (if): a 2 = a 4 = a 6 = and 



a ±^ + 4/1(1 + 2*/) 
a 1 = a, = a 5 = ^—^ , (91) 

F H = ~R 2 [2fi- R 2 (l + 2u)]-aR i . (92) 

(5) Mixed mode, a 2 = a 4 = a 6 = 0, ai = a 3 7^ a 5 . It is always unstable. 

(6) One-dimensional quasipattern, a 2 = a 4 = 0, ai i3j 5 i6 7^ 0. 

(7) Two-dimensional Quasipattern, a\ — • • • qq — R, with R given by, 



_ a ± v /a 2 + 4Ml + 2z/ + 2 7 + 2/3) 

2(l + 2z/ + 2 7 + 2/5) ' ^ 

F Q = -3R 2 [2/2 - i? 2 (l + 2v + 2 7 + 2/5)] - 2ai? 3 . (94) 

The relative stability of the former solutions is as follows: 

(1) Rolls are always unstable at onset. Provided v > 1, (5 > 1/2, and 7 > 1, 
they become stable at a larger value of /x, given by 

(2) Rectangles are unstable at onset and become stable at 

a 2 (l + 7) 

if v > 1 and v + 2(3 > 7 + 1 > 0. 

(3) Squares are also unstable at onset. Ifz/ + 7>l+2/3>0 they become 
stable at 

a = a2{1 + m (97) 
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(4) Hexagons appear in a saddle-node bifurcation at 



a 2 



" = "4(TT2^ (98) 

but when 2(7 + /?) < —(1 + 2v) they are always unstable with respect to 
the dodecagonal quasipattern. If 2(7 + /3) > —(1 + 2v) they are stable 
at the saddle-node but can become unstable as a result of a secondary 
bifurcation. In particular, if v > 1 they can become unstable to rolls at 

and, when 1 + 2v > 2(7 + (3), to the dodecagonal quasipattern for 

" > (i + 2,:VW (100) 



(5) Dodecagonal quasipatterns appear in a saddle-node bifurcation at 

la 2 

/i = ~4(l + 2z/ + 2 7 + 2 / 3) 2 - (101) 

When 2|7 + /3| < 1 + 2u, they are unstable at onset with respect to 
hexagons, but become stable through a secondary bifurcation, at 

a 2 67 + 6^-2^-1 
^Ta + 2,-27-2^- (1 ° 2) 

Therefore, there is hysteresis between hexagons and the dodecagonal 
quasipattern. A similar situation is found in the case of superlattices 
[21]. 

Besides, ifl + 7<z/ + 2/3orl + 2/3<7 + z/ the quasipattern becomes 
again unstable to rectangles or squares, respectively. The values of the 
control parameter for these transitions are 

* iree ~ (1 + 7-^-2/3)2- [L[)6) 



and 



= a 2 (2 + 7 + ^ + ^) 
^ (l + 2/ 9-7 -z/) 2 ' 1 J 
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4-2 Longwave analysis 



The longwave analysis proceeds analogous to that of the previous cases. The 
linearized perturbation equations are given by: 

d t ri = -2qR(hi ■ V)<pi + (n* • V)\ + Ra(r i+2 + r i+4 - n) (105) 
-2R\ - 2uR 2 (r l+2 + r i+A ) - 2 1 R 2 (r l±1 + r l±3 ) - A(3R 2 r lTl , 

dtfa = (hi ■ Vf<Pi + -Ifa • V)n - aRifa + l+2 + i+4 ), (106) 

with 



= a + v /a 2 + 4( / i-^)(l + 2z/ + 2 7 + 2/j) 

2(1 + 2v + 2 7 + 2/3) ' [ ] 



For the amplitude perturbations there are four eigenvalues. Two correspond 
to one dimensional eigenspaces, = —R[2R(1 + 2v + 2 7 + 2(5) — a], with 
eigenvector v T H = [1, 1, 1, 1, 1, l]/v% and (X T2 = — i?[2i?(l + 2i/ - 2 7 - 2(3) - 
a], with = [ — 1, 1, — 1, 1, — 1, l]/>/6- The other four eigenvalues o"t 34 = 
— 2i?[i?(l + 7 - 2/3 - i/) + a] and a T5 6 = -2R[R(1 + 2(3 - 7 - v) + a] are 
associated with two two-dimensional eigenspaces. Four orthonormal vectors 
spanning these spaces are: 

vl = [-1/2, 1/2, 1/2, -1/2, 0, 0], v t Ta = [-1/2, 1/2, 0, 0, 1/2, -1/2], (108) 

corresponding to (Jt 3A and 

vl = [1/2, 1/2, 0, 0, -1/2, -1/2], vl 6 = [0, 0, 1/2, 1/2, -1/2, -1/2] (109) 



to <7r 4i5 . 

The global phases are again stable modes with one-dimensional subspaces 
spanned by u% l = [1, 0, 1, 0, 1, 0]/y/3 and u% 2 = [0, 1, 0, 1, 0, 1]/V3- There are 
four marginal modes. The two phase modes corresponding to translations in 
space can be written as 



ul = [1, 0, -1/2, - V3/2, -1/2, V3/2]/V3, (110) 
ul = [0, 1, y/3/2, -1/2, -V3/2, -1/2J/V3. (Ill) 

An orthonormal base for the phasons, which corresponds to a relative trans- 
lation of the two hexagonal lattices, is given by: 
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< = [1, 0, -1/2, VS/2, -1/2, ->/3/2]/>/3, 
< 2 = [0, 1, -V3/2, -1/2, -l/2]/>/3. 



(112) 
(113) 



As in the case of the octagonal quasipattern, under rotations the phason field 
<p = (ipi, ip 2 ) changes with five times the rotation angle. 

An expansion very similar to that in the decagonal case leads to the longwave 
equations 



d t $ = D x V 2 + (D 2 - £>i) V(V • 0), (114) 
d m = D3VV1 + (£>4 - D 3 )d y (d y ip 1 + d x cp 2 ), (115) 
d m = D 3 V 2 <p 2 + (D 4 - D 3 )d x {d xVl + d y <p 2 ), (116) 

— * 

with = ((fi x , (f) y ). Note that to this order the equation for the phonon modes 
of the dodecagonal pattern is the same as that for an isotropic medium. The 
values of the coefficients are: 



1 q 2 

*=!----■ < 118 » 

4 Vi Ui 

1 q 2 

D, = --!L, (119) 

D. = \-*-£, (120) 
4 w 2 u 2 

with = -i2[a - 2R{2(3 + 2v + 2 7 + 1)], w 2 = -R[a + 2R{2(3 + 2 7 - 2v - 1)], 
«i = 2i?[a + + 7 - 2/3 - 1/)], u 2 = 2R[a + R(l + 2/3 - 7 - 1/)] and i? given 
by Eq. (93). 

In complex form the longwave equations become: 



d t <f> = \(D 1 + D 2 )\ V| 2 + ip 2 - £i)V(V0 + V0), (121) 
= \{D 3 + D 4 )|V|V + ^3 - £> 4 )VV (122) 



One feature that distinguishes the dodecagonal quasipattern from the octag- 
onal and the decagonal quasipattern is the fact that its longwave dynamics 
decouple into pure phonon and pure phason modes. Thus, the eigenvalues are 
simply given by 

a x = -D x Q 2 , a 2 = -D 2 Q 2 , a 3 = -D 3 Q 2 , a A = -D A Q 2 , (123) 
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Fig. 9. Stability diagrams for a = 1 and a) v = 0.7, 7 = 0.3 and (3 = 0.5, b) v = 0.9, 
7 = 0.8 and (3 = 0.4, and c) v = 0.9, 7 = 0.9 and (3 = 0.2. The lower solid line 
corresponds to the saddle-node bifurcation of the quasipattern. In all the figures 
the quasipattern is unstable at onset and becomes stable above the dotted line (cf. 
Eq. (102)). In a) and c) there is a further transition to rectangles and squares, 
respectively, above the upper solid line. The stable region of the quasipattern is 
limited by longwave instabilities corresponding to the phonon (long dashed lines) 
and phason modes (dashed-dotted lines). The diamond and square in c) correspond 
to the simulations in Figs. 10 and 11, respectively. 



and one obtains separate phase and phason instabilities when D 12 or .D34 
change sign, respectively. Typical stability diagrams are shown in Fig. 9. The 
stability limits due to phonon and phason modes are indicated by dashed and 
dashed-dotted lines, respectively. Over wide ranges of parameters the stable 
wavenumber band is limited by phason modes. 

For the phonon-type modes, the eigenvectors correspond to the usual irrota- 
tional and divergence free modes (equivalent to transversal and longitudinal 
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Fig. 10. Instability corresponding to 03 > (diamond in Fig. 9c) for a) t=0, b) 
t=380, and c) t=450. (ji = 5, v = 0.9, 7 = 0.9, (3 = 0.2, q = 0.4, L = 50, 

k c = 16fc m m) • 

waves in an elastic medium.) They satisfy: 

V-0 ; = O, Vx0 t = O. (124) 

For the phason modes, on the other hand, the eigenvectors are: 





Qx 


e lQ - x , 


<p°< = 


Qy 


e iQx 




Qy 






Qx 





(125) 



These modes satisfy the equations: 

d y <p? + d x <p? = 0, d x ip?-d y ip? = 0. (126) 
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Fig. 11. Instability corresponding to 04 > (square in Fig. 9c) for a) t=0, b) 
t=440, c) t=445 and d) t=600. {ji = 20, v = 0.9, 7 = 0.9, = 0.2, q = 0.5, L = 25, 

k c = 16fc m m)- 

In order to study the behavior arising from these instabilities we have simu- 
lated numerically Eqs. (85), (86). We start with a perfect dodecagonal quasi- 
pattern and add a perturbation in the form of (p a ' A (Fig. 10a) or <fi a4 (Fig. 
11a), with Q = (An/L, Att/L). Since the perturbation is along the diagonal, 
the evolution of the system will be quasi one-dimensional. In Fig. 10 we show 
the evolution of the instability corresponding to 03 > 0. The perturbation 
grows until it creates line defects at time t=380 (Fig. 10b) along which four 
of the amplitudes vanish and change their wavenumber, demonstrating that 
the instability is subcritical. In Fig. 10c the final state is shown, after the line 
defects have gone. It corresponds to a slightly distorted quasipattern. 

The instability corresponding to a 4 > is also subcritical, generating line 
defects (Fig. lib). But for the values of the parameters in the simulation the 
square pattern minimizes the Lyapunov functional and it nucleates around 
the line defects. The patches of squares grow (Fig. 11c) until they occupy the 
whole cell, resulting in a slightly distorted square pattern (Fig. lid). 
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5 Defects 



The properties of defects in the quasipatterns differ noticeably in the three 
cases investigated in this paper. Due to the absence of resonance terms in 
(1) (also to higher orders) the octagonal quasipattern can be considered as 
made up of four sub-lattices corresponding to the four basic wavevectors. Each 
of the sub-lattices has its own defects. Thus, taking into account the two 
possible 'charges' of the defects there are 2x4 different defects. They interact 
strongly with defects of the same sub-lattice via the phase. The interaction 
with the defects in the other sub-lattices is only through the variations in the 
magnitude and is much weaker. Overall, we expect the dynamics to be quite 
similar to that of defects in square patterns, which are made of two rather 
than four sublattices. Topologically, two defects of different sub-lattices could 
bind to form a vectorial defect (cf. [26]). However, in our simulations we have 
not observed such kind of defects in the octagonal quasipattern. This seems 
to be consistent with simulations of the Vector Complex Ginzburg-Landau 
Equation, where vectorial defects were never observed in the potential limit 
[27]. 

For the dodecagonal quasipatterns the Ginzburg-Landau equations (85,86) 
have (quadratic) resonance terms coupling the three modes in each of the two 
hexagonal sub-lattices. At the defects therefore two amplitudes in the same 
sub-lattice must vanish leading to penta-hepta defects very similar to those in 
the usual hexagonal patterns. However, while in the core of penta-hepta defects 
of hexagons the pattern corresponds to the roll pattern, it corresponds in the 
dodecagonal case to the one-dimensional quasipattern. As in the hexagonal 
case, each penta-hepta defect carries two independent charges corresponding to 
the two vanishing amplitudes. This leads to a total of (2 x 3) x (2 x 2) x 2/2 = 12 
different defects. While in the octagonal case collisions between defects can 
only annihilate them or leave them unchanged, collisions between penta-hepta 
defects can change their type. An example of such a collision is given by the 
process 



(+,+,0; 0,0,0) + (-,0,+; 0,0,0) - (0, +, +; 0, 0, 0). (127) 

Here +/— in the i th entry stands for a defect with positive /negative charge in 
the mode Ai while indicates that the corresponding amplitude has no defect. 
The semi-colon separates the two hexagonal sub-lattices. While there are a 
number of different such type-changing collisions, a penta-hepta defect in one 
sub-lattice can never change into one in the other sub-lattice, since there is no 
resonance term in (85,86) involving modes of both hexagonal sub-lattices. In 
other words, there is no penta-hepta defect that has one vanishing amplitude 
in one sub-lattice and one vanishing amplitude in the other sub-lattice. 
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Penta-hepta defects within the same sub-lattice interact with each other strongly 
through the phase. The strength of the interaction (attraction vs. repulsion) 
is expected to be related to the sum N of the products of the charges of the 
individual defects [28], 

N = ±5)5l (128) 
i=i 

Here 8j' 2 is the topological charge of the first and second defect, respectively, 
in the mode Aj. For hexagons (n = 3) N can only take on the values -2, -1, 
1, and 2, since for any defect pair there is always a mode that vanishes in 
both of them. Thus, within each sublattice all penta-hepta defects interact 
strongly. Penta-hepta defects of different sub-lattice will interact only weakly 
(through the magnitude). Considering that collisions do not change defect type 
across the sublattices, one may expect the ordering dynamics that leads from 
disordered patterns to ordered quasipatterns to occur in the two sub-lattices 
essentially independently. 

The decagonal case appears to be the most interesting one in terms of the 
defect dynamics. As in the case of defects in hexagon patterns, the (quartic) 
resonance term requires that in a defect two amplitudes vanish. In contrast 
to the hexagonal or dodecagonal case there are, however, two qualitatively 
different defect types. In one of them the vanishing Fourier modes are rotated 
by 27r/5 with respect to each other while in the other they are rotated by 
2 x 2rt/5. While in the former case the core exhibits the quasipattern H 2 , it is 
the quasipattern Hi that appears in the latter case (cf. Fig. 5d,e). Examples of 
the two cases are shown in Fig. 12a,b. Overall, there are (2 x 5) x (2 x 4)/2 = 40 
different defects. 

The same arguments employed in the hexagonal case [28] suggest that the 
interaction between defects is richer in the decagonal quasipattern than in 
hexagons or in the other quasipatterns investigated in the present paper. In 
hexagonal or dodecagonal patterns two penta-hepta defect pairs always share 
a mode that vanishes in both pairs and the charge product iV given by (128) 
never vanishes, implying that the defects always interact strongly through 
the phase. In the decagonal quasipattern, however, the analogously defined 
N (with n — 5) can also take on the value 0, since two pairs of defects need 
not share a mode with vanishing amplitude. The above argument therefore 
suggests that in this case the interaction is very weak. We have confirmed this 
expectation by numerical simulations of the Ginzburg-Landau equations (39). 
While defect pairs that share a vanishing amplitude (N ^ 0) attract or repel 
each other quite strongly, those with N = move so slowly that our prelim- 
inary simulations were not able to identify even whether their interaction is 
attractive or repulsive. 
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Fig. 12. Decagonal quasipattern with a defect of type a) H\ and b) Hi in the center. 
The top panels show a cross section of the five amplitudes in the x-direction. Two 
of them vanish at the core of the defect. Of the remaining three, two are equal, as 
expected from Eqs. (44), (46). In the bottom panels a reconstruction of the physical 
field ij) is shown. At the core of the defects, the one-dimensional quasipatterns H\ 
and H2 can be observed (compare with Figs. 5d,e). 

A separation into weak and strong interaction also arises in the dodecagonal 
quasipattern. There, however, the division between strong and weak is parallel 
to the division into the two sublattices. In the decagonal case, however, the 
resonance term involves all modes and the quasipattern cannot be viewed as 
the combination of two separate sub-lattices. Correspondingly, a defect of one 
type can change into a defect of any other type through collisions with suitable 
other defects. We expect that the exceedingly weak interaction of defect pairs 
with N = may slow down the evolution from disordered or random initial 
conditions towards an ordered quasipattern. 



6 Conclusions 

In this paper we have addressed the stability of various types of quasipat- 
terns with respect to long-wave sideband instabilities. For quasipatterns with 
octagonal, decagonal, and dodecagonal rotational symmetry we have derived 
from the corresponding Ginzburg-Landau equations long-wave equations for 
the two phase and the two phason modes. Their stability analysis yields the 
long-wave stability properties of these patterns. 

For the octagonal and the decagonal quasipatterns the phase and the phason 
modes are coupled and one cannot clearly separate the instabilities in those of 
the phases and those of the phasons. In these cases our numerical simulations 
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suggest that the nonlinear behavior arising from the instabilities is the same 
as that observed in the usual Eckhaus instability. Thus, the instabilities do 
not saturate within the long-wave equations and lead to the creation of defect 
pairs which subsequently annihilate each other yielding a stable quasipattern 
with slightly modified wavevectors. Interestingly, however, in the dodecagonal 
quasipatterns the phase and phason equations decouple and there are param- 
eter regimes in which the quasipatterns first becomes unstable with respect 
to phason modes rather than phase modes. The ensuing dynamics appear to 
be somewhat different than for the usual phase modes. This question has, 
however, not been pursued in detail in this paper. 

The interaction between defects of the decagonal quasipattern can be ex- 
tremely weak for certain combinations of defects. We expect that this will 
have a strong influence on the evolution from disordered initial conditions to 
the regular quasipattern. The investigation of the evolution can most likely 
not be addressed within the Ginzburg-Landau equations (39) since the order- 
ing dynamics will also entail a spreading of the modes in Fourier space along 
the critical circle (cf. the dynamics found in hexagon patterns with rotation 
[29]). This will necessitate the use of suitable equations of the Swift-Hohenberg 
type. 

We gratefully acknowledge discussions with J. Vinals, L. Kramer and M. Sil- 
ber. This research was supported by NASA (NAG3-2113), the Engineering 
Research Program of the Office of Basic Energy Sciences at the Department 
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